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1. Introductory Remarks and Calculation Setup 

We will present the main features of the method and preliminary results of the bag parameter 
calculation for the K meson oscillations at three values of the lattice spacing using the Nf = 2 
dynamical quark configurations produced by the ETM collaboration. 

ETMC dynamical configurations have been produced with the tree-level Symmanzik improved 
action in the gauge sector while the dynamical quarks have been regularized by employing the 
twisted mass (tm) formalism [|l|]. It has been demonstrated that with the condition of maximal twist 
this formalism provides automatic 0(a) -improved physical quantities 

In the so called physical basis the fermion lattice action concerning the sea sector is written 

Ssea = ^ V/(x) ( /V - iys T3 Wcr + pL.ea ) , ( 1 ■ 1 ) 

X 

with Wcr = — f V* +Mcr{r = I); Y = {u d)^ is a doublet of degenerate light sea quarks while 
l^sea = diag(/i„ /!(/). We should also note that the tm formalism offers a simpler renormalisation 
pattern with comparison to the standard Wilson regularization. This is true for some important 
physical quantities calculated on the lattice, as for example the pseudoscalar decay constant and 
the chiral condensate. 

It has been shown that the use of the tm regularization can simplify the renormalization pattern 
properties of the four-fermion operators which enter in the calculation of certain phenomenolog- 
ically important weak matrix elements such as Bk [H, H, Q]. In order to achieve both 0{a) im- 
provement and a continuum-like renormalization pattern in the evaluation of Bk we introduce the 
valence quarks with Osterwalder-Seiler lattice action and allow for replica of the down (d, d') and 
strange {s, s') flavours viz. 

S„ai=a'^Y^ qf{x)iyy -iy^rfWcr + jlAqfix) , -r^ = = rj/ = = 1 . (1.2) 

X f=d,d'.s,s' 

The valence sector action above is written (unlike eq. ( |1.1| )) in the so called physical quark basis 
with the field qf representing just one individual flavour. While the four fermion operator relevant 



for Bk (see eq. (2.1)) is chosen to contain all the four valence flavours in eq. ( 1 .2 ), the interpolating 
fields for the external (anti)Kaon states are made up of a tm-quark pair (dy^s, with — = r^) and a 
OS-quark pair (d'y^s', with r^/ = r/). This mixed action setup with maximally twisted Wilson-like 
quarks has been studied in detail in Ref. allows for an easy matching of sea and valence quark 
masses and leads to unitarity violations that vanish as as the continuum limit is approached. In 
the present case, however, the quark mass matching is incomplete because we are neglecting the 
sea strange quark (i.e. partially quanched computation), thereby inducing some (possibly small) 
0(a°) systematic error. We notice that the proposed method for obtaining automatic 0{a) improved 
results has already been tested successfully in the calculation of Bk with fully quenched quarks |^]. 

In Table | we give the simulation details concerning the mass values of the sea and the valence 
quarks for each value of the gauge coupling for the calculation presented in this work. The smallest 
sea quark mass corresponds to a pion of about 270 MeV for the case of j8 = 3.90. For j8 = 4.05 
the lightest pion weights 300 MeV while for j8 =3.80 the lowest pion mass is around 400 MeV. 
The highest sea quark mass for the three values of the lattice spacing is about half the strange quark 
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mass. For the inversions in the valence sector we have made use of the stochastic method (one- 
end trick of ref. [0]) in order to increase the statistical information. Propagator sources have been 
located at randomly chosen timeslices. For more details on the dynamical configurations and the 
stochastic method application see Refs ^. 
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Table 1: Simulation details 



2. The K-meson bag parameter 



We recall that in our mixed action setup all the physical quantities are evaluated with no 0(a) 
discretization effects (see Ref. [^) and moreover the four fermion operator relevant for Bk, which 
reads 



phys-basis 
bare 



[{qsYiiqd){qs'Yiiqd>) + {qsrn75qd){qs'7tir5qd')] + [d^ d'] . (2.1) 



is multiplicatively renormalizable. This can be easily understood by noting that in the (unphysical) 
tm basis, where the Wilson term enters the valence action in the standard way (with no //j -twist) 
and the operator renormalization properties are the same of the standard Wilson fermionic action, 
the operator ( 2A ) takes the form 

tm-basis 
bare 



V^A^ +A^V^ 



[iXsY^Xd)iXs' 7^1 YsXd' ) + iXsYti Y5Xd)iX'^' Y^iXd' )] + [d^ d'] , (2.2) 



Here Xf = exp ''^'^^''^qf and Xf = ^/exp ^y^^l'^J = d,d' ,s,s' are the tm basis valence quark fields. 
The operator (2.2) is known to be protected from mixing under renormalisation due to CPS sym- 



metry [10]. In summary we have ("R" stands for "renormalized") 

phys-basis 



V^V^+AaA 



phys-basis 
R 



WA+AV 



V^V^ +AuA 



bare 



^VA+AV 



tm-basis 

bare 
(2.3) 



where the name of the renormalization constant is chosen so as to be consistent with the notation 
used in the standard Wilson fermion literature. 

In order to estimate the Sx-parameter we calculate a three-point correlation function where a 
four-fermion operator is free to move in lattice time t whereas two "K-meson walls" consisting of 



3 



K-bag parameter from tmQCD 



S. Simula 



noisy sources are imposed at fixed time separation tR — ti = T/2. The ti value changes randomly 
from configuration to configuration. In our simulations we consider the time reversed case too and 
we average them properly. The plateau signal is taken for ti <^ t <^ tR. We extract from the 
ratio: 

In our analysis all correlation functions satisfy the condition apLi = a/Xsea while the valence strange- 
Uke quark mass values are given in Table |l[ An important remark is in order: the mixed regu- 
larization set-up that we have used leads at finite lattice spacing to different values for the decay 
constant and the pseudoscalar masses of the two K-mesons employed in the calculation. We find 
that the discretisation effects are negligible for the decay constant while happen to be significant in 
the case of the pseudoscalar mass. For this reason we normalize the four fermion matrix element by 
dividing with (8/3)/m^^ m'jl f^^ f^. Moreover, as expected, the cutoff effects diminish drastically 
towards the CL. So this kind of systematic error is well under control. 

The fits to the light quark mass behaviour are performed using the SU (2) Partially Quenched 



Chiral Perturbation Theory formula of refs [ ]11| , |12| ]. In our case the fit ansatz is: 



B{^h) = B^i^n) [l +Z;(At/O^At/ - 3^2i"'ln(^)] +D{^H)a^ (2.5) 

where jXh denotes the quark mass values around the strange quark (see Table |l|). Thus, the fit 
procedure consists of a combined fit of chiral and continuum extrapolation. We find that the cutoff 
effects on our data are well described by a /X/ -independent (but /X/, -dependent) 0{a^) term. 

Two methods of analysis have been followed. The first method relies on using the information 
for the physical mass values of the up/down and strange quarks in the continuum limit, as they have 



been estimated in a recent ETMC computation [13]. Note that the implementation of this method 



requires the knowlegde of the quark mass renormalization constant []14p. The second method con- 
sists of employing the pseudoscalar masses instead of the quark masses. In this case we choose 
a set of three values of reference pseudoscalar masses made out of two strange-like quarks, M/,/,; 
keeping each of them fixed we perform the chiral fits in terms of the light pseudoscalar mass. In 
the end of the procedure we estimate Bk via an interpolation at the physical point defined by the 
formula Mj^. = 2M\ — Both methods give compatible final results within less than one standard 
deviation. 

In Figure [I]^a) the quality of the plateau is shown for j8 = 3.90, for three values of the light 
quark mass and for one typical value of /i/,; in Figure |I|(b) we present an example of a combined 
chiral plus continuum fit (three value of the lattice spacing) for B^^{l,h) versus the light pseu- 
doscalar mass squared in units of ro; the value of M^h is in the vicinity of the physical one. 

The two point renormalisation constants for the axial and vector current have been calculated 



using the RI-MOM method [15]. We recall that the physical axial current made up of OS quarks is 



normalized by Za while the one consisting of tm quarks is normalized by Zy [14]. The RI-MOM 



method has also been employed for the calculation of the renormalisation constant of the four- 



fermion operator [16]. In Figure 0(a) we show the behaviour of the renormalisation constant as a 



function of the momentum squared in lattice units {apY for j3 = 3.90 at the valence chiral Umit 
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Figure 1: (a) The quality of the plateau for three values of the light quark mass for j3 = 3.90; (b) Com- 
bined chiral and continuum fit for B^^{l,h) versus (roM//)^. The empty black circle gives the value at the 
continuum limit for the case of (rOM/,/,) = 1 .50; 



and for afJ-sea = 0.0040. Discretization effects of 0{a^) have been evaluated at one loop [ [T7| ] and 
subtracted from the relevant correlation functions. Thus, the leading discretization effects on our 
RI-MOM determination of the renormalization constant are of 0{g^a^,g^a^). We show three types 
of results; two of them correspond to two estimates of the subtracted perturbative contributions. 
The amount of the subtraction depends on the choice of the value for the gauge coupling. We have 
considered two cases for the gauge coupling, the naive (go) and the boosted one (gh). We also 
show the result for the Zy^_|_^v/(RGI) without considering any perturbative subtractions (indicated 
as "uncorrected" in the figure). In the right panel of the same figure we illustrate the absence of 
mixing with "wrong chirality" operators; in fact the mixing coefficients are vanishing. 





P a 

(b) 



Figure 2: (a) RI/MOM computation of the multiplicative renormalization factor Zy^+AV (RGI) at fi — 3.90; 
(b) Mixing coefficients (k ~ I, - ■ ■ ,4) with other four-fermion operators with "wrong chirality". 



Our preliminary result for Bk in the RGI scheme in the continuum limit is 

B|gi = 0.73(3) (3) 
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The first error includes the uncertainty coming from the correlators and from the fit procedure 
(chiral plus continuum) while the second one is due to the uncertainties in the calculation of the 
renormalisation constants. We are currently attempting to reduce the latter uncertainty. 

3. The K-bag meson parameter beyond the SM 



Interactions beyond the SM including supersymmetry furnish new diagrams in the calculation 
of the AS' = 2 process. Their effect expressed in the OPE expansion is to enrich the set of the 
local operators to be considered in the low energy regime, see e.g. JT^]. Therefore one has to 
calculate on the lattice the matrix elements of five parity even four-fermion operators namely 0\ = 
Ovv+AA, O2 = Oss+pp, O3 = Ofj,04 = Oss-pp and O5 = Ow^aa E ||, 

It is well known that the renormalisation pattern of the parity-even four-fermion operators 
becomes complicated because of mixings as soon as the regularization breaks the chiral symmetry; 
this is certainly the case of Wilson fermions. However using the proposal of Ref. [Q] this problem 
is bypassed; due to the axial rotation mapping of the parity-even to parity-odd operators in the tm 



basis the renormalisation pattern becomes continuum-like [16]. It is worth mentioning that, as in 
the case of the SM four-fermion operator, the lattice estimates of the matrix elements of O2, • • • j O5 
are automatically O (a) -improved. 

First results regarding the signal quality for the case of j3 = 3.90 are given in Figure p|. We 



depict the plateaux for the B3 bag parameter (left panel) and for the quantity 



{K\Oi\K) 



(right 



panel). Both figures refer to the same value of the light quark mass for three different choices of 
the strange-like quark mass. Computation at the other two values of the lattice spacing as well as a 
full determination of the renormalisation constant matrix is still in progress. 
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Figure 3: (a) and (b): the quality of the signal for the quantities B3 and Rt, respectively at three values of the 
strange-like quark mass using a/X/ = 0.0040 at j3 = 3.90. 
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